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^ . Abstract 

c . 

^ , If the neutrino mass is non-zero, as hinted by several experiments, then i?-parity-violating supersym- 

metric Yukawa couplings can drive a heavy neutrino decay into lighter states. The heavy neutrino 
' may either decay radiatively into a lighter neutrino, or it may decay into three light neutrinos 

fT^ , through a Z-mediated penguin. For a given mass of the decaying neutrino, we calculate its lifetime 

for the various modes, each mode requiring certain pairs of i?-parity-violating couplings be non-zero. 
' We then check whether the calculated lifetimes fall in zones allowed or excluded by cosmological 

^ , requirements. For the latter case, we derive stringent new constraints on the corresponding products 

' of _R-parity-violating couplings for given values of the decaying neutrino mass. 

^ ! PACS number(s): 13.35.Hb, 12.60. Jv, ll.30.Fs 

o 
o 

o , 

"Th ! 1 Introduction 



Oh! 



The Super-Kamiokande collaboration has recently provided compelling evidence in support of neutrino 
oscillations as an explanation of the atmospheric anomaly The observed solar neutrino deficit 
(— I I and the LSND accelerator experiment results are also indicative of oscillations. These results achieve 

especial significance since oscillations require the neutrinos to be massive. A massive neutrino signals 
physics beyond the standard model and has far-reaching implications in particle physics, astrophysics, 
and cosmology B. 



X 



H I In this work, we consider i?-parity-violating (^) supersymmetry (defined later) in the context of 

neutrino decays ^. Such interactions violate lepton number and in the presence of appropriate 
couplings of this type a heavier neutrino of one flavour can decay to a lighter one of a different flavour in 
association with the emission of a photon. Alternatively, the heavier neutrino can decay invisibly into 
three lighter neutrinos through one loop graphs involving ^ couplings. 

Cosmological and astrophysical requirements forbid certain regions in the neutrino mass and lifetime 
plane. They originate, for example, (i) from the precise black-body nature of the microwave background 
radiation, (ii) from the tight requirements of consistency of the predicted primordial nucleosynthesis 
with observation, etc., which are discussed later. 

We calculate, for a heavy neutrino of a given mass, the decay lifetime induced by ft couplings. 
If the lifetime so obtained falls in the forbidden region then couplings of the chosen strength are not 
allowed for the neutrino mass used. In this way we can establish new constraints on some lepton number 
violating couplings applicable for specific neutrino masses. On the other hand, if one assumes that these 
^ couplings are at their existing limits, then upper bounds can be set on the decaying neutrino mass. 

Before we move to the next section, a short introduction to i?-parity and ^ supersymmetry is in 
order. 'i?-parity' in supersymmetry refers to a discrete symmetry which follows from the conservation of 
lepton-number (L) and baryon- number (B). For the standard model particles and their superpartners. 
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it is defined as i? = (^_iY3B+l+2S) ^ where S is the intrinsic spin of the field. R is +1 for all standard 
model particles and —1 for all superparticles. 

The most general ft superpotential is given by, 

Wif = \\,kUL,El + KjkUQM + ^^l.^W'jDl + ^Ji,UHu, (1) 

where fc — 1,2,3 are quark and lepton generation indices; Li and Qi are iS'C/(2)-doublet lepton 
and quark superfields respectively; Ef, J7f, Df are S'C/(2)-singlet charged lepton, up- and down-type 
quark superfields respectively; iJ„ is the Higgs superfield responsible for the generation of up-type 
quark masses; Xijk and X'^^j. are L-violating while A^^j, are i3-violating Yukawa couplings. B- and L- 
conservation are not ensured by gauge invariance and hence there is a priori no reason to set these 
couplings to zero. A^j^ is antisymmetric under the interchange of the first two generation indices, while 
^'ijk antisymmetric under the interchange of the last two. Thus there could be 27 A', 9 each of A and 
A" couplings and 3 /ii parameters. We assume that the generation indices correspond to the flavour basis 
of fermions. We also note at this point that the B-violating couplings A-'^-^ and the bilinear couplings 
Hi are not of any relevance to our present analysis. We shall deal only with the L-violating trilinear 
couplings and, for the sake of simplicity, we assume that they are real. 

Stringent constraints on individual L-violating couplings have been placed from the consideration 
of neutrinoless double beta decay, ^'e-Majorana mass, charged-current universality, e — fi — r universality, 
deep-inelastic scattering, atomic parity violation, r decays, D and K decays, Z decays, etc. Product 
couplings (two at a time), on the other hand, have been constrained by considering — e conversion, 
jj, 67, b — > S7, B decays into two charged leptons, — Ks and Bq — Bq {q = d, s) mass differences, 
etc. (For a collection of all these limits, see |^). 

In the following section we elaborate on the neutrino mass spectra which are preferred by the 
experimental data. In section 3 we discuss the radiative neutrino decay modes driven by the A'- and 
A-type ^ couplings while in section 4 we turn to the invisible three neutrino decay channel. In the 
next section we check if and how cosmological requirements on these decay lifetimes constrain the ^ 
couplings. We end in section 6 with our conclusions. 



2 Neutrino mass spectrum 

In this section we elaborate on the neutrino masses that we will be using to examine possible neutrino 
decays driven by ^ interactions. We are motivated in making these choices by the available evidence 
for massive neutrinos. 

The upper limits on neutrino masses obtained from laboratory experiments are the following. The 
masses of ly^ and Vt are constrained to be rUi,^ < 0.17 MeV and m^^ < 18.2 MeV As regards Ve, the 
upper limit from tritium beta-decay end-point measurements is 2.5 eV |l^, while if it is of Majorana 
nature then from the absence of neutrinoless double beta decay one has the constraint^ < 0.2 eV 

Before we discuss the constraints that emerge from the oscillation data, a few remarks are in order. 
The ^ couplings have been defined in terms of flavour eigenstates. On the other hand, any discussion 
of decay must refer to a neutrino in its mass eigenstate. For simplicity of presentation and the ease of 
illustration, in the following we do not distinguish between the flavour and the mass eigenstates. Indeed, 
the indications for neutrino oscillation imply that this is not the actual situation. In fact, the mismatch 
between the mass and flavour bases is an essential ingredient of neutrino oscillations and even maximal 
mixing is sometimes preferred (e.g., as required by the atmospheric anomaly). It is straightforward to 
incorporate the effect of this mixing in our results. One must multiply the decay rates for Vi —^ Vi>, 
presented later, by appropriate factors determined by the probability of the Vi {lyi') being present in the 
parent (daughter) neutrino. 

In much of the analysis we choose the decaying neutrino to be and for the purpose of illustration 
take its mass to range from 45 eV to 100 keV. This is consistent with the laboratory bound on m^^ 
We consider its decay to either i/^ or 1/^. We also consider the possibility of i^f^ decaying to i/g. There are 

^Strictly speaking, this bound applies to tlie (ee) element of the Majorana mass matrix written in the flavor basis. 
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models in the literature proposing inverted mass hierarchy scenarios m^^ < m^^ < m^^. But, as far as 
bounds on the .^-couplings are concerned, these provide no new constraints since the rate for i^i — > 
decay flavour indices) in this scenario is the same as that for — > Vi if the masses m^. and m^., 
are interchanged. 

Now let us recollect what is known about the magnitude of neutrino mass splittings from the 
data on neutrino oscillations. There are two independent recent results - namely, the atmospheric 
neutrino anomaly |Q] and the solar neutrino problem - which can be conveniently explained within 
the framework of neutrino oscillations. The LSND experiment has also claimed positive evidence of 
i^fj. — i^e (and also v^ — v^^j oscillations However, this result is still awaiting independent confirmation. 

Oscillations require a non-zero mass splitting between neutrinos of different flavour. More precisely, 
the experimental data determine Am^ = m\ — rn^, where r7ii,2 are the masses of two neutrinos. Oscilla- 
tions also determine the unitary transformation which relates the mass basis of neutrinos to the flavour 
basis. As mentioned earlier, we are not directly concerned with this aspect and do not discuss it any 
furthei|5 

The information from experimental results on neutrino mass splittings are as follows: 

1. The atmospheric neutrino data can be explained in terms of oscillations of to either the Vr or 
a sterile state v^. The sterile neutrino does not couple to any of the standard model particles and 
passes undetected through the experimental set up. The necessary mass splitting has been found 
to be Am2 - f 0"^ - cV^ = Aatmos- 

2. The solar neutrino problem can be explained in terms of vacuum neutrino oscillations of the v^. to 
either the v^,t- or the Vg with Am^ ~ 10~^° eV^. The data admits an alternative solution in terms 
of MSW resonant neutrino conversion if Am^ ~ 10~^ — 10~^ eV'^. In the latter case the sign of 
the mass splitting is fixed and the v,, must be the lighter state of the two. We indicate the mass 
splitting required for a solution to the solar neutrino problem by Agoiai - 

3. The LSND experiment indicates — Vf, oscillation with Am^ ~ 1 cV^ = Alsnd- 

As is seen from the above, if both the atmospheric neutrino anomaly and the solar neutrino deficit 
are sought to be explained by invoking oscillations between three sequential neutrinos only, then the 
mass splittings among them must be small. On the other hand, if one admits a fourth (sterile) neutrino 
in the model, then it is possible to arrange for larger mass splittings. With a fouth neutrino, one can 
have two categories of mass hierarchies. In category 'A', there are three closely spaced almost degenerate 
states explaining the solar and atmospheric data and the fourth one is separated by a larger gap A. 
The fourth neutrino could either be the sterile or v.,- since the Ve {i^p.) must be closely spaced with some 
other state to satisfy the solar (atmospheric) constraint. In category 'B', there are two pairs of almost 
degenerate neutrinos separated by the larger gap (details later). Taking the LSND experiment into 
consideration, 'A' is incompatible with the data, while 'B' can be comfortably accommodated |l3| . 

Two cases may arise if we work in 'A'. The heaviest neutrino could be (dominantly) sterile. Since it 
does not have any direct ft couplings, its decay via mixing with other states will be strongly suppressed, 
and we cannot set any bound on the ^ couplings. Therefore, this case is not interesting for our purpose 
and we will not consider it in our analysis. On the contrary, the heaviest state could be (dominantly) 
the z/t, decaying via ^ interactions. Since our analysis concerns, after all, the decay of a heavy active 
neutrino into at least one lighter active state, this latter case does not yield any different result from 
what we obtain if we work in 'B'. Hence, to simplify matters, we adopt the framework 'B' for our 
subsequent discussions. 

In 'B', which is a 'two-pair' scenario, all four neutrinos participate in the solution to the solar 
and atmospheric neutrino anomalies. In this case there are two allowed forms of the mass spectrum 
consisting of two closely spaced pair of states (spacings Aatmos and Agoiar) with a larger separation 
(A) between the two pairs as shown in Fig. 1. Inclusion of the LSND result fixes the mass splitting 
A to Alsnd fl^- It needs to be mentioned that the oscillation results (excepting for the MSW case) 

^ In this work we motivate possible neutrino mass spectra from the experimental results and do not attempt to reproduce 
the mass matrix structures that would fit the data. There is a vast and growing literature in this subject related to ft 
supersymmetry pd]. 
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determine only the magnitude of the spacing between the levels. Therefore, in addition to the two mass 
spectra shown in Fig. 1, other possibilities obtained from these by exchanging the members within any 
of the two closely spaced pairs or by interchanging the relative ordering of the two pairs gives rise to 
spectra which are equally acceptable. As alluded to earlier, our results below are not affected by such 
inversion of hierarchies. 



A. 



)lar 



A 



^solar 



(a) 



(b) 



Figure 1: The two allowed forms of the neutrino mass spectra from atmospheric and solar results. Note that 
the spacings Aatmos and Agoiar are only indicative and not to scale. If the LSND results are also included then 

A = Alsnd. 

In our subsequent work, pending confirmation from another independent experiment, we do not 
impose the LSND constraint on the mass spectra. Once we abandon the LSND results, both 'A' and 'B' 
are acceptable by the remaining data. In any case, as noted before, except for omitting the uninteresting 
case of the decaying heavy sterile neutrino in the framework 'A', we do not lose any generality by 
sticking to the option 'B|J. Thus, we assume that there are two pairs of states with Arn^ determined 
by the atmospheric and solar neutrino results but we leave the splitting between these pairs arbitrary, 
demanding only that consistency with the laboratory limits be maintained. This permits us to consider 
mass splittings as large as 170 keV (but not more), which is the laboratory upper limit on m^^, and check 
the cosmological implications if such heavy neutrinos decay into lighter states via the |? interactions. 
It should be noted that, in this framework, the laboratory upper limit of 18.2 MeV on m^^ can never 
be reached, as Fig. la suggests that m^^ can at most be a little above m„ . The generation indices for 
the decaying heavy neutrino and the lightest neutrino produced as a decay product are chosen as the 
combinations (3,1), (3,2), and (2,1) respectively. We will neglect the mass of the product neutrino. 

3 Radiative decay of a heavy neutrino 

A heavy neutrino, which from now on we denote by , can decay radiatively into a lighter neutrino 
(anti- neutrino) , denoted by v (i/), through AL = (2) penguin diagrams containing ^ couplings. We 
discuss these two cases separately. 

3.1 The decay ^ u + 'j 

The AL — decay of h'^ to v and a real photon occurs through penguin diagrams (see Fig. 2). Although 

^We make a remark in passing that Big-Bang Nucleosynthesis (BBN) data constrain the mixing between an active and 
a sterile neutrino pj]. For example, the constraint on the additional number of neutrino species ATV^ < 0.2 15 disfavours 
the large angle Ufj^-i/g mixing solution to the atmospheric neutrino problem. A 'safer' estimate AA^^ < 1 |l6) , on the other 
hand, can accommodate an arbitrarily strong mixing between an active and a sterile state. For the present analysis we do 
not indulge ourselves further along these lines. 
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Figure 2: Penguin diagrams corresponding to the AL = radiative decay of i^^. 



the right-handed neutrinos do not couple in the ]^ superpotential, a mass insertion on the heavy neutrino 
Hue (which requires to have Dirac-type mass) enables one to write down the following effective 
Hamiltonian for this decay: 

1. ^ .... - 



where 



Here 



M - e KjkK'jkQdNc-^-T-^ {Ilk - Iij) ■ 



(3) 



(4) 



with r^j. = [maJrhY'. In the above expression, and in all subsequent discussions, we assume a common 
mass m for whichever scalar is exchanged. Qd is the charge of the down-type quark inside the loop. 
A^c = 3 is the colour factor. 



3.2 The decay z/^ ^ z/ -|- 7 

This AL — 2 decay of a, v^j to v and a real photon occurs through penguin diagrams (see Fig. 3). 
Contrary to the previous mode, for this decay to take place the decaying heavy neutrino does not 
require to have a Dirac mass. The effective Hamiltonian for this decay is given by 



AL=2 



where 



Here 



A2 = e X^jk\'k]QdN, 
1 



{X2k —llj)- 



L23 



+ln(rdj) 



(5) 
(6) 
(7) 



As is seen from Fig. 3 this mode involves the mixing between left- and right-type d-squarks. This 
has been approximated by Am^ (LR) = {A — ^jl tan /3) md ~ fhmd in the above expression. 
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(b) 





(d) 




Figure 3: Penguin diagrams corresponding to the AL = 2 radiative decay of v^^. 

It follows from Eqs. (||) and ^ that for the special case when j = k, the radiative decay amplitude 
vanishes identically. A look at Figs. 2 and 3 reveals that for each diagram in which a photon is attached 
to a particle (fermion or scalar) of generation j there exists a similar diagram in which the photon 
couples to the corresponding anti-particle of generation k - hence the cancellation when the generations 
match. 

In the above we have considered processes which are triggered by non- vanishing A' couplings. These 
processes could also be driven by A-type couplings. The transition can be readily performed in Figs. 2 
and 3 by replacing the internal down-type quarks with charged leptons and the down-type squarks with 
sleptons of the corresponding generations. Indeed, Nc would be unity in this situation. 



4 The invisible decay of into three Ught neutrinos 

The other kinematically allowed decay mode of is into three neutrinos/anti-neutrinos. Such a decay 
proceeds via a set of Z- mediated penguin diagrams (see Fig. 4). The choice of A' (or A) couphngs 
determines the flavour of one final state neutrino. In fact, for the aforenoted case oi j — k leading 
to vanishing radiative decays, such three-body invisible channels constitute the only decay modes of 
Vfj. Since the Z boson, unlike the photon, couples to fermions (sfermions) and their conjugates with 
different strengths, the cancellation observed for a radiative decay process does not occur for a Z- 
mediated penguin. As we will see later, the decay widths in these three-body modes are much smaller 
than the two-body radiative decay widths. In the following, we shall consider only the AL = penguins. 
The AL = 2 diagrams suffer a suppression from left-right sfermion mixing, and we will not consider 
them in our subsequent discussions as their numerical impact is insignificant. It is also possible to drive 
this decay through box graphs involving the products of A couplings, but their contribution to the same 
order as that of the penguins, we have checked, is insignificant. 

For simplicity, we present the expressions for the AL = penguin for those cases in which the 
internal fermion and the sfermion are of the same type (i.e., j = k) and the heavy neutrino mass is 
negligible in comparison with the internal fermion masses. In such a situation the one- loop effective 



vertex can be written as 



:(5aL7pL'L, where, 



rd, Td^ ln(rd, ) 



(8) 
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5 Cosmological constraints 



Cosmology sets tight constraints on allowed neutrino masses and associated lifetimes. For a neutrino 
of a certain mass, any decay mode driven by ^ interactions has a lifetime determined by the couplings 
and the exchanged particle masses. It needs to be checked whether such decays are consistent with the 
cosmological requirements (discussed below) or whether one obtains new bounds on the ^ couplings. 
We discuss the different modes in the following subsections. In order to permit easy comparison with the 
existing bounds on such couplings, we choose, as in the literature, all exchanged supersymmetric scalar 
particle (squark or slepton) masses to be 100 GeV in the discussion below. The bounds corresponding 
to any other sfcrmion mass can be readily obtained by using the formulae given in this paper (although 
just by scaling arguments one can derive these numbers to a very good approximation). 

Though we are examining neutrino decays, the case of a cosmologically stable neutrino plays an 
important role. If the lifetime of the neutrino is larger than the age of the universe (^o) [ pT|, t hen 



it is cosmologically stable. They must either be rather light (the Cowsik-McClelland bound [18|) or 
comparatively heavy (the Lee- Weinberg bound jl^). The allowed values for a stable neutrino mass 
are ||] 

< 92h^no eV, or > ^ \ GeV (9) 

Here, JIq is the ratio of the present energy density to the critical energy density and h is the normalized 
Hubble expansion rate, h is defined as Hq = lOO/i km s~^ Mpc~^, where Hq is the Hubble expansion 
rate. Present data indicate that ft,^f2o is smaller than unity. Since we are considering neutrinos of mass 
100 keV or less, the second inequality in Eq. (^) is not of relevance for us. For the sake of discussion 
we take the first limit to be 45 eV; its exact value depends on the precise magnitude of the combination 
h'^rio- Hence, a neutrino of mass 45 eV or more must necessarily be cosmologically unstable. 

5.1 The — > u(u) + 7 mode 

The final state photons in radiative decays of neutrinos are subject to constraints from the cosmic 
microwave background radiation (CMBR) spectrum, primordial nucleosynthesis, etc. For the AL = 2 
radiative process ^ P + 7, the neutrino decay lifetime is 

"''^^^ ^ K jeV)]3 [A'i']^ 

where /i lies in the range 4.8 x 10^^ — 1.4 x 10^^ depending on the masses of the quarks and squarks 
exchanged in the loops. [A' A'] denotes the product of the relevant ^ couplings. When the A couplings 
are involved, /i lies between 1.0 x 10^^ — 3.6 x 10^^. 

For the sake of numerics, we choose a reference value of rrii,^ = 45 eV. Since perturbativity requires 
all ^ couplings be smaller than unity, the lifetime of the neutrino from Eq. ( p^ is larger than 5.3 x lO^'^s, 
comparable to tree ^ 3 x lO^^s (the recombination epoch). Now t^, ~ tree for a 45 eV neutrino is not 
admissible from cosmological requirements on the following grounds. If to > T,y > tree, then the photons 
produced in the radiative decay undergo redshift and contribute to the diffuse photon background since 
after recombination no charged particles are available to scatter them. This sets a bound ||20|| : 

{^)"'<z.,o'w,'r" {^y\ (11) 

A 45 eV neutrino cannot satisfy this constraint if its lifetime is as estimated above. 

The remaining possibility then is t^, > to. These photons would simply get superimposed on the 
diffuse photon background. The experimental data, in this situation, yield the bound: 

The above equation can be translated to upper bounds on the ^ couplings. In Table 0, we have collected 
these new upper bounds emerging from the AL = 2 radiative decay for a 45 eV decaying neutrino. Side 
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Combinations 


Existing bounds 


New bounds 


X' X' 


1 


X 


10-3 


1 


X lO-'i 


X' X' 
^321 ^112 


7 


X 


10-3 


6 


X 10-3 


X' X' 

^323^132 


1 


X 


10-1 


1 


X 10-5 


X' X' 

'^331'^113 


3 


X 


10-3 


1 


X 10-'' 


X' X' 

^332*^123 


7 


X 


10-3 


1 


X 10-5 


X' X' 

'^313'^231 


9 


X 


10-3 


1 


X lO-'' 


X' X' 
'\321'^212 


2 


X 


10-2 


6 


X 10-3 


X' X' 
^^323^^232 


1 


X 


10-1 


1 


X 10-5 


X' X' 
'\331'^213 


9 


X 


10-3 


1 


X 10-'* 


X' X' 
^332^^223 


1 


X 


10-2 


1 


X 10-5 


X' X' 
'^213^^131 


1 


X 


10-3 


1 


X 10-'^ 


X' X' 

^223^^132 


2 


X 


10-2 


1 


X 10-5 


X' X' 

'^231^113 


4 


X 


10-3 


1 


X 10-4 


X' X' 

^232^123 


1 


X 


10-2 


1 


X 10-5 


A323A132 


3 X 10-3 


5 X 10-^ 


A323A232 


3 X 10-3 


5 X 10-5 


A232A123 


2 X 10-3 


5 X 10-5 



Table 1: New upper bounds on products of different A' and A couplings for a iieavy neutrino mass of 45 eV and 
for AL = 2 radiative decay. The existing bounds are obtained by multiplying the upper bounds on the individual 
couplings Q. 



by side, we have displayed the existing constraints on these couplings. We have presented only those 
cases for which our new bounds are more stringent than the existing ones and have rounded off these 
numbers to the first significant digit. 

Alternatively, one can turn the argument around and find the maximum mass of the decaying heavy 
neutrino that is acceptable if the R couplings achieve their present experimental upper limits. We have 
presented these bounds in Table ^ using the most dominant mode — the AL = 2 radiative decay. 

Now consider a neutrino of mass of 100 keV (this is near the maximum limit of 170 keV, as we argued 
in section 2). Since ^-couplings must be smaller than unity, one obtains, from Eq. (|l0|), t^,^ 4.8 x lO'' s. 
But dumping of extra photons by such massive neutrinos decaying at a time close to the thermalization 
epoch and much after nucleosynthesis is not cosmologically acceptable |2^, ^ , the relevant bound from 
the black-body nature of the CMBR being: 

For the AL — radiative decay — > 1^ + 7, the amplitude is proportional to the Dirac mass of 
the decaying neutrino. So the behaviour for the lifetime is as follows: 

"''('^ ^ KjeV)]5 [A'A']2 ^2, (14) 

where /2 hes between 2.1 x 10''° — 1.7 x 10^5^ por the A couplings, /2 hes in the range 4.2 x lO** — 
9.9 x 10^5, Por a 45 eV decaying neutrino, Eq. (^4|) implies a lifetime larger than the age of the universe 
even after setting the relevant ^ couplings at their experimental upper limits. Thus, i?-parity violating 
couplings cannot be constrained from this process and if the only non-vanishing ^ couplings turn out 
to be those that drive this AL — decay, then a neutrino of mass 45 eV, in order to be cosmologically 
acceptable, must have other interactions resulting in a faster decay. It is easy to convince oneself from 
the scaling argument that for a decaying neutrino mass as high as 100 keV, the lifetime is either larger 
than the age of the universe or, if less, certainly larger than tree, the recombination epoch. In the 
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Combinations Existing bounds Corresponding maximum 

in cV 



X' X' 

313 131 


1 X 10"^ 


15 


X' X' 

321'^112 


7 X 10-3 


40 


X' X' 

^323^132 


1 X 10-1 


0.4 


X' X' 

331 ^113 


3 X 10-3 


8 


X' X' 

^332^^123 


7 X 10-3 


2 


X' X' 
^313^^231 


9 X 10-3 


5 


X' X' 
^321^^212 


2 X 10-2 


26 


X' X' 

^323-^232 


1 X 10-1 


0.4 


X' X' 

^331^^213 


9 X 10-3 


5 


X' X' 
^332^^223 


1 X 10-2 


1 


X' X' 
^213^^131 


1 X 10-3 


15 


X' X' 
^223^132 


2 X 10-2 


1 


X' X' 

^231^113 


4 X 10-3 


8 


X' X' 

^232^^123 


1 X 10-2 


1 


A323A132 


3 X 10-3 


6 


A323A232 


3 X 10-3 


6 


A232A123 


2 X 10-3 


7 



Table 2: Maximum allowed values of the heavy neutrino mass decaying in the AL — 2 radiative mode, based on the 
assumption that the A' and A couplings achieve their existing upper limits. These existing limits are as in column 2 
of Table 1. Recall that in Table 1 we have set new upper limits on the relevant product couplings assuming a mass 
of 45 e\/ for the decaying neutrino. 



latter case, the neutrino mass and radiative decay lifetime must satisfy relation (^_l|), which is clearly 
not possible. Therefore, a 100 keV neutrino decaying via ^ couplings alone is inadmissible. 

5.2 The ^ u + u + 1? mode 

Usually the radiative decay mode dominates over the other loop decay modes and consequently the latter 
are of any numerical relevance only if the former is absent. The radiative decay amplitude involves the 
combinations ^ijkK'kj fo'" ^ v + j mode (see Eq. (||)) and X'ijkK'jk the ^ v + j process 

(see Eq. (H)). In the special case j — k, the amplitude vanishes identically in each case, as can be easily 
seen from Eqs. (|) and (|). The same situation prevails for the A couplings. Therefore, for bounds on 
the KjjK'jj ^ijj-^i'jj) combinations we must turn to other processes. It is in these situations that 
the present mode becomes of relevance. 

For the AL — Z^-mediated z^^^ — + I'l'D process (the two neutrinos coupled to the Z-vertex must 
have the same flavour) the lifetime is given by 

"-^^^ ^ K jeV)]MA-Af 

Here lies in the range 10^2 — 10^2 Pqj. ^jjg \ couplings, the range is Iff''' — 10^^. It is seen from 
Eq. ( p^ that for a 45 eV (or even 100 keV) neutrino, the lifetime for this process is always larger than 
the age of the universe, to. Therefore, the |?-couplings cannot be constrained from this process for such 
neutrinos. 
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6 Conclusions and Discussions 



We conclude by highlighting the salient features of our analysis. The main thrust is to constrain from 
cosmological considerations those ^ couplings which trigger a heavy neutrino decay either radiatively 
into a lighter neutrino or into three lighter neutrinos. Such attempts were also made in the past ^. 
Our analysis, though similar in nature, differs with the previous ones in the neutrino mass spectrum 
used and complements them by including a relevant piece of the Lagrangian overlooked in the existing 
literature. 

There are two issues that one needs to address. One is the origin of the neutrino mass and the 
other is the interaction that drives its decay. We have considered a four-neutrino (three standard and 
one sterile) framework which permits mass-splitting (or for that matter even the mass of the decaying 
neutrino) as large as order 100 keV or so, remaining consistent with the oscillation data and the recent 
result of the tritium beta-decay experiment. The authors of ref. |^ have restricted themselves to a 
three-neutrino framework and considered neutrino masses not more than order 1 eV. In fact, as we 
argued in section 2, a three-neutrino scenario, given the present constraints from oscillation and tritium 
beta-decay, cannot permit a neutrino to have a mass more than order 1 eV. In this situation it is not 
possible to improve on the existing constraints on the ^ couplings. Obviously the four-neutrino model 
provides some breathing space allowing at least one neutrino to have a significantly larger mass decaying 
(with more phase space) to the lighter species. The resulting constraints, depending on the mass of the 
decaying neutrino, are indeed tighter. In the analysis of refs. ||, ^, the neutrino Majorana masses are 
generated only by ^ interactions and the same couplings eventually lead to their decay. In our case, on 
the contrary, the neutrinos receive masses (either Dirac- or Majorana-type) from some other origin. In 
fact, we have put these masses in by hand, holding the ^ interactions accountable only for their decay. 
If instead one assumes that the ^ interactions responsible for neutrino decay also generate the mass of 
the decaying neutrino then the limits we obtained will be altered. This can be seen as follows. The 
combinations ^ij^K'tj ^^Sut have been constrained by our present analysis also contribute to the off- 
diagonal (ii') terms of the left-handed (neutrino) Majorana mass matrix, while the combinations ^'ij^^'n^j 
contribute to the diagonal (ii) term. Assuming that m^^ in Eq. (|lO|) is generated in this manner and, 
to make a simple estimate, choosing the off-diagonal {ii') element of the matrix about a tenth of this 
diagonal term, we can derive an upper bound on A^^j.A^,^^ demanding consistency with the cosmological 
constraint in Eq. (p^). In this way, the strongest limit that we obtain is A323A232 < 10~^. This indicates 
that if ^ interactions alone are responsible both for the generation of neutrino mass as well as for its 
decay, then cosmological constraints imply that the mass so generated would be comparable to 45 eV. 

The other point we would like to stress upon is that from the X'^-^^LiQjD'^ term in the superpotential 
we get two kinds of squark-mediated processes, one involving the scalar from Qj and the other the one 
from D'^ . They correspond to two different pieces in the Lagrangian. For the radiative decay there is an 
exact cancellation between these two contributions for j = k. This point has been overlooked in ^. 
It is in this situation that the Z-mediated penguins provide the only decay mode of the heavy neutrino. 
Similar arguments can be advanced for the A couplings. 

It should be noted that the product couplings which drive the ~^ Vi decay are also responsible for 
a similar process involving their SU(2) partners, namely, the Z^/ — > li decay (where I is a charged lepton). 
We have checked that the constraints we have derived on the product couplings in the context of neutrino 
decays are consistent with the bounds on the same combinations derived from the non-observation of 
the charged lepton decays . 

The vvv mode which we have discussed arises from a flavor- violating Z-v^-v coupling 

generated by penguin diagrams involving i?-parity violating couplings. Such a non-diagonal Z vertex 
can also be produced at tree level due to (a) neutrino-neutralino mixing arising from the misalignment 
between the /i^ and the sneutrino vacuum expectation values, or (b) through mixing between active and 
sterile neutrinos. The magnitude of these mixings depend on details of the model which are beyond 
the scope of the present analysis. However, one can readily see that the lifetime of a 45 eV neutrino 
decaying through such an interaction is approximately T7, (m^/m,y^)^ sin^ (/) (here sint/) parametrizes the 
suppression at the off-diagonal Z vertex), which is much larger than the age of the universe. As a result, 
such decays can never compete with the radiative decay modes which yield the bounds which we have 
presented. 



11 



It should be noted that neutrino osciUations will not affect the results shown in Tables 1 and 2. 
Since the decay times that we find are larger than the age of the universe while the typical oscillation 
times, as required by solar and atmospheric data, are at most a few minutes, there is little scope of 
repopulation of the vh states from the daughter neutrinos through oscillation. However, a very large 
number of oscillations will have taken place before the decay so that the decaying states will have an 
averaged effect, the probability of an initial vh remaining so being cos^ 6, where is a mixing angle. If 
to start with an equal number of neutrinos of all types exist then this averaging will leave this balance 
unaltered. 

Before closing, a few comments on the existing constraints from baryogenesis ||2^ are in order. 
The requirement that GUT-scale baryogenesis does not get washed out imposes A" ^ 10"^ for generic 
indices, although it has been argued that such bounds are model dependent and can be evaded [^ . 
In our analysis we have assumed all A" to be zero - this also solves the proton decay problem. In the 
absence of the A" couplings, the A' or A couplings alone cannot wash out the initial baryon asymmetry, 
unless there are B violating but [B — L) conserving sphaleron-induced non-perturbative interactions. 
Since the latter interactions conserve {^B — Li) for each lepton generation, the assumption that not all 
lepton numbers are simultaneously violated suffices to preserve the initial baryon asymmetry. Keeping 
this mind, whenever in the present work we have considered a product of two A' (or A) thereby violating 
two lepton numbers, we have made an implicit assumption that the third lepton number is conserved. 



Acknowledgements 

GB acknowledges hospitality of the Max-Planck Institute in Heidelberg, and CERN Theory Division, 
Geneva, where some parts of the work were done. SR acknowledges support from the Council of Scientific 
and Industrial Research, India. AR has been supported in part by the Council of Scientific and Industrial 
Research and the Department of Science and Technology, India. 



References 

[1] Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998). 

[2] B.T. Cleveland et al, Nucl. Phys. B (Proc. Suppl.) 38, 47 (1995); Kamiokande Collaboration, Y. 
Fukuda et al, Phys. Rev. Lett. 77, 1683 (1996); GALLEX Collaboration, W. Hampel et al, Phys. 
Lett. B 388, 384 (1996); SAGE Collaboration, J.N. Abdurashitov et al, Phys. Rev. Lett. 77, 4708 
(1996); J.N. Bahcall and M.H. Pinsonneault, Rev. Mod. Phys. 67, 781 (1995). 

[3] LSND Collaboration, C. Athanassopoulos et al, Phys. Rev. Lett. 75, 2650 (1995); ihid. 77, 3082 



(1996); nucl-ex/9706006 



[4] See, for example, Rabindra N. Mohapatra and Palash B. Pal, Massive Neutrinos in Physics and 
Astrophysics, 2nd Ed., World Scientific, Singapore (1998). 

[5] K. Enqvist, A. Masiero, and A. Riotto, Nucl. Phys. B 373, 95 (1992). 

[6] E. Roulet and D. Tomniasini, Phys. Lett. B 256, 218 (1991). 

[7] G. Farrar and P. Fayet, Phys. Lett. B 76, 575 (1978); S. Weinberg, Phys. Rev. D 26, 287 (1982); N. 
Sakai and T. Yanagida, Nucl. Phys. B 197, 533 (1982); C. Aulakh and R. Mohapatra, Phys. Lett. 
B 119, 136 (1982). 



For recent reviews, see G. Bhattacharyya, Nucl. Phys. B (Proc. Suppl.) 52A, 83 (1997); 



ph/9709395, Invited talk given at Workshop on Physics Beyond the Standard Model: Beyond 



the Desert: Accelerator and Nonaccelerator Approaches, Tegernsee, Germany, 8-14 Jun 1997; H. 



Dreiner, hep-ph/9707435 , published in 'Perspectives on Supersymmetry', Ed. by G.L. Kane, World 



Scientific, Singapore; R. Barbier et al. Report of the Group on i?-parity violation, |hcp-ph/9810232 



12 



[9] Particle Data Group, C. Case et at, Eur. Phys. J. C 3, 1 (1998). 

[10] The TROITSK Collaboration, V.M. Lobashev et ai, Phys. Lett. B 460, 227 (1999); The Mainz 
Collaboration, C. Weinheimer et ai, Phys. Lett. B 460, 219 (1999). 

[11] L. Baudis et ai, Phys. Rev. Lett. 83, 41 (1999). 

[12] S. Rakshit, G. Bhattacharyya and A. Raychaudhuri, Phys. Rev. D 59, 091701 (1999); R. Adhik ari 
and G. Omanovic, Phys. Rev. D 59, 073003 (1999); A. Abada and M. Losada, |hep-ph/9908352| ; G. 
Bhattacharyya, H.V. Klapdor-Kleingrothaus and H. Pas, Phys. Lett. B 463, 77 (1999); M. Drees, S. 
Pakvasa, X. Tata and T. ter Veldhuis, Phys. Rev. D 57, 5335 (1998); E.J. Chun, S.K. Kang, C.W. 
Kim and U.W. Lee, Nucl. Phys. B 544, 89 (1999); O.C.W. Kong, Mod. Phys. Lett. A 14, 903 (1999); 
A.S. Joshipura and S. Vempati, Phys. Rev. D 60, 111303 (1999); A. Datta, B. Mukhopadhyaya and 
S. Roy, Phys. Rev. D 61, 055006 (2000); B. Mukhopadhyaya, S. Roy and F. Vissani, Phys. Lett. B 
443, 191 (1998). 

[13] S. Goswami, Ph ys. Rev. D 55, 29 31 (1997); S.M. Bilenky, C. Giunti and W. Grimus, Eur. Phys. J. 
C 1 247 (1998); |hep-ph/971 13111 C. Giunti, |hep-ph/ 990 7485 . 



[14] A.D. Dolgov, |liep-ph/0006103 . 



[15] D. Tytler, J.M. O'Meara, N. Suzuki and D. Lubin, |astro-ph/0001318 



[16] E. Lisi, S. Sarkar and F.L. Villante, Phys. Rev. D 59, 123520 (1999). 

[17] According to the Big Bang model of cosmology tg is estimated tobell.5 + l±1.5 Gyr [^. 

[18] R. Cowsik and J. McClelland, Phys. Rev. Lett. 29, 669 (1972); S. Gerstein and Ya.B. Zeldovich, 
Pisma Zh. Eksp. Teor. Fiz 4, 174 (1972). 

[19] B.W. Lee and S. Weinberg, Phys. Rev. Lett. 39, 165 (1977); P. Hut, Phys. Lett. B 69, 85 (1977); 
K. Sato and M. Kobayashi, Prog. Theor. Phys. 58, 1775 (1977); M.I. Vysotsky, A.D. Dolgov and 
Ya.B. Zeldovich, Pisma Zh. Eksp. Teor. Fiz. 26, 200 (1977). 

[20] E.W. Kolb and M.S. Turner, The Early Universe, Addison- Wesley, Redwood City (1990). 

[21] See Fig. 2 of G. Gelmini and E. Roulet, Rept. Prog. Phys. 58, 1207 (1995). 



[22] G. Gelmini, |hep-ph/9904369| ; S. Sarkar, Rept. Prog. Phys. 59, 1493 (1996). 

[23] K. Huitu, J. Maalampi, M. Raidal and A. Santamaria, Phys. Lett. B 430, 355 (1998); M. Chaichian 
and K. Huitu, Phys. Lett. B 384, 157 (1996). 

[24] A. Bouquet and P. Salati, Nucl. Phys. B 284, 557 (1987); B.A. CampbeU, S. Davidson, J. Ellis and 
K. Olive, Phys. Lett. B 256, 457 (1991); W. Fishier, G.F. Giudice, R.G. Leigh and S. Paban, Phys. 
Lett. B 258, 45 (1991). 

[25] H. Dreiner and G.G. Ross, Nucl. Phys. B 410, 188 (1993). 



13 



